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Abstract 

The regularization of the CMB power spectrum is an important issue of cosmology. Most 
approaches assume that there is no need to regularize the power spectrum, while Parker 
advocated the new regularization approach for the power spectrum in 2007 [1]: the adiabatic 
regularization, which was originally developed for particle creation in curved spacetime (see 
the review [2]). This thesis focuses on this issue, especially concerning adiabatic subtraction 
terms. The subtraction terms for minimally coupled slow-roll inflation are well known (see 
e.g. [1]). We extend the view to more generic inflation models, and derive the subtraction 
terms for k-inflation models. Via the method of Urakawa-Starobinsky [3], we consider the 
time development of the subtraction term at late times. We also consider the non-minimally 
coupled case, and show that the adiabatic regularization is independent of the frame: Jordan 
or Einstein frames. 
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Part I 

Introduction 


The inflationary paradigm [4-7] is a fascinating approach to solving various cosmological 
problems. Although inflation can describe the solutions naturally, we do not know how it 
occurs, or even whether or not it really occurred. 

To examine inflation, we have a tool to see the early universe: the cosmic microwave 
background (CMB) radiation. The CMB consists of photons which became free at early 
times (the last scattering surface) and now fill the observable universe. We can observe the 
temperature and the polarization, and we know that the temperature of these photons from 
each direction is almost identical. 

However, tiny fluctuations in the temperature and the polarization have a large amount 
of information of the early universe because it is considered that they reflect the primordial 
perturbations of inflaton or gravitation. If we can analyze the observation data properly, it 
means that we can see the universe at the last scattering surface. Indeed, the CMB is the 
oldest light that we can observe. To know the details of the universe before the last scattering 
surface directly, we are forced to detect other signals such as gravitational wave signatures. 

We can infer from the CMB that the early universe went through an inflationary era [8]. 
However, the precision is still not enough to determine which inflation model was realized, 
and many inflation models have been advocated from various theoretical motivations. To 
compare theoretical inflation models with the observed CMB, the power spectrum which will 
be shortly explained in section (1.2) is used. Therefore the method used to estimate the 
theoretical power spectrum is an important matter of cosmology. 

Although the amplitude of the perturbation which is the integrand of the power spectrum 
has divergences (before regularization), usually it has been assumed that regularization is not 
necessary 1 . However, since the paper of Parker [1] quite a lot of discussion has arisen con¬ 
cerning this issue. He applied the adiabatic regularization method, which is a regularization 
method used in quantum field theory in curved spacetime, to the power spectrum. 

The purpose of this thesis is to review the adiabatic regularization for the power spectrum 
and to expand it to more general cases. In the rest ot this part, Part 1, we review cosmological 
perturbation theory and we also review adiabatic regularization and k-inflation in the first 
half of Part 2. Then we expand our interest to non-minimally coupled k-inflation, and derive 
the adiabatic subtraction terms for the model in the last part of Part 2. In Part 3, we discuss 
the issue of the pliysicality Einstein and Jordan frames that are related by a conformal 
transformation, where of relevance to this thesis we show the conformal equivalence of the 
adiabatic subtraction terms in each frame. In Part 4 we conclude and mention future work. 

We include detailed calculations in Appendices: In Appendix A, we calculate Hubble flow 
functions in detail. In Appendices B and C, the time development of the adiabatic subtraction 
terms is discussed in some specific models. Finally, in Appendix D, we try to express the 
subtraction terms in the Jordan frame and discuss the generality of the derived subtraction 
terms. 

1.1 Cosmological perturbation theory 

1.1.1 The Friedmann equations 

Observationally, the universe is nearly homogeneous and isotropic and hence perturbation 
theory is usually considered around such a background. In this thesis, we go up to linear 
order in perturbations, although in principle higher order could be considered. 

1 Recently, a review of the problem [9] has even claimed that there is no need to regularize the power spectrum 
from the standpoint of observations (although interacting theories would still require a regularization according 
to [9], see also [10]). 


4 



The background metric 

The homogeneous and isotropic background universe is written by the Friedmann-Lemaitre- 
Robertson-Walker (FLRW) metric. 


ds 2 


dt 2 
dt 2 


a 2 {t)dx 2 


a 2 (t) 


dr 2 

1 - Kr 2 


+ r 2 {d9 2 + sin 2 dd(/) 2 ) 


( 1 . 1 . 1 ) 


where K is spatial curvature. a{t) is a scale factor which denotes the size of the universe at 
specific time. In this thesis, only the spatially flat (K = 0) case will be considered. 

To simplify the calculations, conformal time r) defined as 


V = 



dt' 
ait') 


( 1 . 1 . 2 ) 


is also used. The metric is written in terms of i] as follows 


ds 2 = a 2 {if) [drj 2 — dx 2 ] 


(1.1.3) 


Both proper time t and conformal time i) are used in this thesis. The derivatives is written 

by 

fit) = j t f{t ), f\rf = ±-f{r,) . (1.1.4) 

However, we use Eq. (1.1.3) as the background metric unless otherwise specified. 


Metric perturbations 

Then let us consider the perturbations. The metric perturbations has three modes. 1) Scalar 
perturbation cf>(r), x), ip{g, x), B{i], x), E{i], x). 2) Vector perurbation F(r],x)i,S{r],x)i. 3) 
Tensor perturbation h(r],x)ij. 


$ 9 oo{x) = 2a 2 cj) 

6goi{x) = a 2 {B ri +Si) 

Sgij{x) = a 2 {2if)5ij + 2E,ij+Fij-{-Fjj + hij) 

If we pick all scalar modes only, the line element is written as 

ds 2 = a 2 {g) (1 + 2 (jf dr] 2 + 2 B,i drjdx 1 — j(l — 2ip) Sij — 2 E,ij j dx l dxt 


(1.1.5) 

( 1 . 1 . 6 ) 
(1.1.7) 


( 1 . 1 . 8 ) 


The above metric will eventually be used to derive an equation of motion for the inflation 
perturbations, known as the Muhkanov-Sasaki equation. 


Hydrodynamical perturbations 

Classically, the universe is homogeneous and isotropic, and matter can be approximated by 
a perfect fluid. 

The energy-momentum tensor of perfect fluid is given by following form. 

= {E + P)u,_ l u l , — Pg^v (1.1.9) 

where E is the energy density of the matter and P is the pressure, u^ix) is a four vector 
which satisfy the classical part u ctl = (a, 0, 0,0) and = 1. 

In perturbation theory, E , P and u M also have the spatial-dependent perturbation part. 
Let us write the classical part like E c . e.g.) E{x) = E c {rf + 8 E{tj,x). The perturbation of 
Ui can be decomposed scalar part and vector part, so write them explicitly. 

Sui = 5u\\ ,j +5u±i ( 1 . 1 . 10 ) 

where d l 5u±i = 0 and 5u\\ is the scalar part. 
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( 1 . 1 . 11 ) 


The classical part of the energy-momentum tensor becomes 

T co = E c, T° cl = 0 , Tl ] = -P c 5\ j 

and the linear perturbation part (only scalar modes) are given by 

ST° 0 = 6E , 

ST o = 1 (E c + P c )Su lln , (1.1.12) 

ST\ = —SPS i j . 

These equations will be used later on to derive the Muhkanov-Sasaki equation after we have 
first defined the cosmological background. 


The Friedmann equations 

Using these metric and energy-momentum tensor, we can get the Friedmann equation from 
the Einstein equation R^ = T^J 2 . 

The classical part is 

= l = ~l {Ec + 3Pc) {LL13) 

where H = -. The continuity equation is also derived 

E c + 3H(E C + P c ) = 0 . (1.1.14) 


Or in terms of conformal time, 

a 2 E c 
3 


n 2 = 


a 

a 


= - — (E C + 3P C ) 


(1.1.15) 


E' c + 3 H{E C + P C ) = 0 


(1.1.16) 


where R = — = aH. 

a 

The equation of scalar perturbation is more complicated because the metric perturbation 
contains non-physical degrees of freedom. To deal with this problem, we fix the gauge or 
compose gauge invariant quantities as discussed in the following section. 


1.1.2 Gauge invariance 

Gauge fixing can remove the non-physical degrees of freedom in the perturbed metric Eq. 
(1.1.8). However, the variables in the gauge are not always physical quantity, i.e., they may 
be not gauge invariant. 

In this thesis the conformal-Newtonian gauge will be considered. 


The conformal-Newtonian gauge 

The conformal-Newtonian (longitudinal) gauge is defined by 

ds 2 = a 2 (r]) [(1 + 2(f>) dg 2 — (1 — 2^) dx 2 ] . (1.1.17) 

i.e., B = E = 0. 

In this gauge, the perturbations of the metric and inflaton are gauge invariant. Thus we 
can rewrite metric Eq. (1.1.17) or equations in conformal-Newtonian gauge in terms of gauge 
invariant variables. Using the gauge invariant perturbations $ and T, the metric becomes 

ds 2 = a 2 (r]) [(1 + 24>) dg 2 — (1 — 24/) dx 2 ] . (1.1.18) 

$ corresponds to the Newtonian potential. There seems to be two degrees of freedom for the 
scalar mode. Actually, 4> and 4^ are related via the Einstein equation and not independent. 
In minimal coupling case, these two are identical $ = T. 

For other perturbations, the overlines will be used to express gauge invariant variables, 
e.g., Sip(r],x) -> 6<p(r], x). 

2 In this thesis V8nG = Mp 1 1 is normalized so that Mpi = 1. 
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The comoving gauge 

The comoving gauge which is another gauge fixing way is also considered in some references. 
We will not use it in this thesis, but briefly explain it. 

The comoving gauge is defined by the condition 

0 = 0,0^11=0. (1.1.19) 

The scalar part of the metric becomes 

ds 2 = a 2 (if) [drj 2 — e 2n 5ijdx l dx J ] (1.1.20) 

where 1Z is the comoving curvature perturbation which is a gauge-invariant quantity com¬ 
posed by scalar perturbations in the metric and the inflaton. 

This gauge is particularly well suited to performing calculations using the ADM formalism 
[11] in cosmology, e.g., see [12, 13]. 


1.1.3 Mukhanov—Sasaki equation 

After gauge fixing, the first order Einstein equations can be solved. We can calculate it in 
conformal-Newtonian gauge, and then replace the perturbations by the gauge invariant ones. 

In the minimal coupling case, we can combine the Einstein equations and get the Mukhanov- 
Sasaki equation describing the behavior of the comoving curvature perturbation 7Z. 

vl + {elk 2 - —) u fc = 0 (1.1.21) 


where 


v k = zTZ k , 


i 2 V , I+P 

cM ’ 


dP 
dE ' 


( 1 . 1 . 22 ) 


The comoving curvature perturbation 1Z = HI— [13] is a gauge invariant scalar which 
is conserved and is equal to the curvature perturbation on uniform-density hypersurfaces £ 
on superhorizon scales. 

Finally, we should mention that c s is the ’’sound speed,” the propagation speed of the 
perturbation. In models which have canonical kinetic terms we have c s = 1. 


1.2 The power spectrum 


1.2.1 Quantization 


In this section, the definition and the physical property of the power spectrum of the CMB 
temperature is introduced. 

First, we consider how to quantize the Mukhanov-Sasaki variable v(x) to v(x). At first, 
define the canonical momentum t r, 


_ dC 

n= di 

where C is the Lagrangian. We then assume that the quantized variable 
Mukhanov-Sasaki equation. 

We now require the canonical commutation relation on v and tt: 


( 1 . 2 . 1 ) 


v also obeys 


[v{r],x),v{r], y)\ = [n(r], x), Tr(r), y)} =0 

[v(r],x),7r(n,y)\ = iS(x — y) 

Let us consider the Fourier components 


v{i],x) 


d 3 k 

(2rr)i 


A^v k (r])t 


A l v k{i 7) f 


— ik•: 


with the normalization 

v* k v' k - v* k 'v k = -i . 


( 1 . 2 . 2 ) 

(1.2.3) 


(1.2.4) 


(1.2.5) 
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Then Ar and At have the canonical commutation relation 


\A%, A P ] = [A\, 4 ] = 0 , [A it 4 ,] = S(k - k') . 


k ’ k f 


( 1 . 2 . 6 ) 


1 1 

z 2 (27r) 3 
f dk sin (k\x — y\) k 3 2 


d 3 k\v k \ 2 e in < s ~® 

\n k \ s 


k k\x — y\ 27r 2 


If we set x = y , 


<*’<■>■*» - / T■ 

On the other hand, from Eq. (1.2.7): 

(n 2 (ri,x)) = ^ J dk fc 2 ^(fc) ■ 


(1.2.7) 


The power spectrum of two point correlation function for scalar perturbation V-jz(k) is 
defined by 

(ntfl % ,) = 5(k + k')V n {k) ■ 

Using v k = zTZk and Eqs. (1.2.4)—(1.2.6): 

x)TZ(r], y)) = \{v(v,x)v(v,y)) 


( 1 . 2 . 8 ) 


(1.2.9) 


( 1 . 2 . 10 ) 


:.V K (k) = \U k \ 2 (1.2.11) 

The integrand of Eq. (1.2.9) A^(fc) = z{k) is called the dimensionless power spec¬ 
trum. In this thesis, we will sometimes call A 2 the power spectrum. However, we make it 
clear which we indicate from the context when the distinction is needed. 


1.2.2 Behavior of the power spectrum 

Given that the comoving curvature perturbation 1Z is conserved after the Hubble horizon 
crossing, the power spectrum also time-independent after crossing. This is the reason why 
we are especially interested in large scales. The modes in large scales exit early in the history 
of the universe, and they keep more information of primordial perturbations than the modes 
in small scales. 

Is it known that the dimensionless power spectrum of the scalar perturbation in de Sitter 
spacetime can be exactly solved and it does not depend on k. Moreover, the observations 
support the nearly k -independent power spectra[8]. If the dimensionless power spectrum 
can be expressed (or approximated) by power series, the integral Eq. (1.2.9), which is the 
variance of the perturbation, has the divergences on UV scales. Thus we need to regularize 
the power spectrum in theoretical arguments. 

In next part, we will discuss the adiabatic regularization which is a way to regularize the 
power spectrum. We will see the behavior the adiabatic subtraction term in the latter half 
of Part 2. 
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Part II 


Regularization of the power spectrum 

To obtain the finite expression of power spectrum, adiabatic regularization [1, 14], which is 
one of regularization schemes of quantum field theory in curved spacetime, is considered. In 
this regularization, we perform WKB-like expansion of the solution of perturbations, and 
subtract the isolated divergence part. Then we regard the regularized power spectrum as the 
corresponding one to the observable power spectrum. 

In this chapter, we review the adiabatic regularization method. And then we apply it to 
the power spectrum for k-inflation model. 


2.1 Adiabatic regularization 

Let us consider a scalar field (j>{x) obeying the following equation of motion 3 . 

<t>"{x) + f{ri)(j)'{x) + [g(r})did % + h(rj)] <j>(x) = 0 (2.1.1) 

where /(??), g(rj) and h(g) are functions of time and these can be derived from the given 
metric and Lagrangian 4 . 

In this case, the coefficient functions from the metric depends on only time because we 
are considering up to first order of perturbations. For the power spectra, we will consider the 
equation of motion of the perturbation part only. The background metric is Eq. (1.1.18), 
which has the both dependence of time and spatial coordinates. However, the dependence 
on the spatial coordinates appears as the perturbations. Therefore the spatial-dependence 
of the metric appears the second or higher order perturbation parts, and we can use the Eq. 
(2.1.1) in our discussion without loss of generality. 

Consider the Fourier transformation: 



The following normalization of <j> is also imposed. 

« - = ~i (2-1.3) 

If we define Xk (v) as 

</>&) = Xk(v) exp (ik-x-^ J f{r])dr^j , (2.1.4) 

Xk(jl) obeys the following equation. 

Xk(v) + tt 2 {r))Xk{il) = 0 (2.1.5) 

where 

= ~g(v)k 2 + h{rj) - */ 2 (??) - \f'(v) ( 2 - 1 . 6 ) 

This is an equation of motion of a harmonic oscillator if we regard £l(rj) as the frequency. 
Indeed, in the case of a massive scalar field in the Minkowski spacetime, 


f(v) = 0, g(il) = -1, h{ri) = ml (2.1.7) 

and Eq. (2.1.5) becomes 

Xkiv) + {k 2 + ml) Xk(v) = 0 • ( 2 . 1 . 8 ) 

3 We formulate the adiabatic regularization method using the conformal time for some convenience. The 
original formulation [1, 14] uses the proper time, and the both subtraction terms are equivalent. We can obtain 
the subtraction term in terms of proper time by using either the differential equation by proper time or by the 
variable transformation of the subtraction term derived by Eq. (2.1.1). 

4 We assume that the potential is not steep. 
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(2.1.9) 


k 2 + m ^ = ujf, is the square of the frequency and the solution is 


Xk(v) oc exp(-*w fe ? 7 ) . 

In a curved spacetime, the square of effective frequency ui 2 (r]) = —g(jf)k 2 +h(rf) depends 
on time and the coefficient function f(rf) from the Cliristoffel symbol in covariant derivative 
exists. Therefore we perform a WKB-like (adiabatic) expansion for the solution. 
Introducing a fictitious parameter T in the metric. 


9nv{x) gn„(x/T) 


( 2 . 1 . 10 ) 


The order of differentiation of the metric becomes the adiabatic order of the term. We now 
require the adiabatic condition, i.e., In lowest adiabatic order, Xk should have the form 


1 


Xk 


\/2 uj k (ri) 


exp —i 


u}k(r/)drj‘ 


( 2 . 1 . 11 ) 


Under the adiabatic condition, we can get uniquely the adiabatic expansion of the solution 
Xt, and the solution is schematically written by 


Xk(il) = 


V^W 2n {r,) 


exp l —i 


W 2n (i 1 ')dr 1 ') +0(T 


n—2n—2\ 


( 2 . 1 . 12 ) 


up to 2nth adiabatic order. W 2n is obtained later by the effective frequency of the equation 
of motion each order. 

</>£ become the complete orthogonal set of solutions, and A g and At becomes annihilation 
and creation operators after requiring the canonical commutation relations 


[A* = [ 4 > A U = °< A U = 5 $ *'), 

and an adiabatic vacuum state |0) is constructed. 

Ag|0) = 0 (for all k) 


(2.1.13) 


(2.1.14) 


Neglect the higher order terms in Eq. (2.1.12). i.e., Consider the 2nth adiabatic order 
solution. 

exp (-if W 2n (r)')dri'') (2.1.15) 


(2n) / \ 

Xk W = 


y/2 W 2n ( V ) 

X { k n) ( 7 ?) obeys the following equation. 


x!L {2n) + 


WJ n - WJ n 


2 (W 2 A 
\n \ vv 2n J 


x{ 2n) = 0 


Rewrite Eq. (2.1.5) to 


Xk + 


W 2 - W* 


5 (\V~^ 


Xk = 


W 2 - W 5 


I (w~^ 


- n- 


Xk 


(2.1.16) 


(2.1.17) 


In 2nth adiabatic order, we can choose W so that the left hand side vanishes by using Eq. 
(2.1.16). 


0 = 


W 2 - W 


2 (W ^ 


- ir 


Xk 


(2.1.18) 


i.e., The right hand side is higher adiabatic order than 2nth order. We can get from Eq. 
(2.1.18) 


w 2 = n 2 + w? 


s (w-i'j 


(2.1.19) 


where W in arbitrary order. 

Note that U 2 of Eq. (2.1.6) can be decomposed to the terms of 0th adiabatic order and 
2nd order because f(rj) is 1st adiabatic order. Therefore we can write U 2 in terms of the 0th 
adiabatic order part w 2 and the 2nd order part cr = — q/ 2 (»?) — tfid)- 

n 2 =ujl + c r (2.1.20) 
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W is decomposed as following by using the Eq. (2.1.19). 


W = W {0) + W {1) + W {2) + ■ ■ ■ (2.1.21) 

W (0) = w fc 

W (1) = 0 

w(2) = \“k 5 (W a ) +l u k 1(J 

The unregularized (“bare”) two point function of <f> becomes 


(01 </>(£) 0 ( 5 ) 10 ) = 


and take limit x —> x. 

(O|0 2 (x)|O) 


1 

(27r) 3 
1 

(2t r) 3 
1 


d\d\' (j)^{x)(j)*^, (x)S(k — k') 


d 3 kc/)^x)(j) Ux) 


(2tt) 3 exp (ff(r))drj) 


1 


1 


(2tt) 3 exp (ff(r])drj) 

1 1 


d 3 fce ^(5-5) Xfe ( a . )x *(^) 

j d 3 k\xk{x)\ 2 
dk k 2 \xk(x)\ 2 


To regularize this, we need W 1 up to second order. 

w= 1 


tr(°) + ww + 


2tt 2 exp (/ f(r])drf) 

,o second order. 

= ly-To) + tfo-T 1 ) -|_ h- 


( 2 . 1 . 22 ) 


(2.1.23) 


(2.1.24) 


!E- 1(0) = u;- 1 , W~ 1W = 0, W~ m = -w" 2 lfo (2) (2.1.25) 

Now the divergent part is uniquely isolated by the adiabatic expansion. Because only 
W~ 1 ^ oc k~ x and oc k~ n (n > 3) contain the divergent parts, we subtract these two 

terms and define the physical two point function (minimal subtraction scheme) 


(O|0 (Z)|0)p h y s 


1 1 
27t 2 exp (f f{j])drj) 


dk k 2 


IXfe(^)| 2 


I jlfo-T 0 ) + 1T _1 ( 2 )| 

(2.1.26) 


From Eq. (2.1.25), it is found that the Oth order subtraction term W _1 i°l corresponds to the 
zero-point energy and its time dependence is due to the expansion of the universe and/or the 
non-canonical kinetic term. The second order term does not exist in Minkowski spacetime. 

We can use this method to regularize the power spectrum. The calculation for usual 
slow-roll inflation model has done, and it is shown that the subtraction terms for the scalar 
perturbation and the tensor perturbation become sufficiently negligible after the horizon 
crossing [3] (see also [15]). 

While the adiabatic subtraction terms for slow-roll inflation model is studied very well, 
for other inflation models is not. Hence we generalize it for k-inflation model which has more 
generic Lagrangian than slow-roll model. 


2.2 Minimally coupled k-inflation 

2.2.1 k-inflation 

k-inflation (kinetically driven inflation) is a inflation model advocated by [16, 17]. In this 
model, the inflaton have non-canonical kinetic terms and they make it possible that the 
universe expands exponentially. We will review briefly k-inflation model in this section. 

The Lagrangian of inflaton in k-inflation model given by 

= P(tp, X ) = -V(ip) + K(y)X + L(v)X 2 + ■■■ (2.2.1) 
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where ip is inflaton and X = ^d^ipd^p. However, the explicit form of the Lagrangian will 
not be specified except Appendix B, and the general form P(ip, A') will be used. 

The action in minimally coupled case is 

S=IJ<Pxy/=9[-R + A)] . (2.2.2) 

The energy-momentum tensor under the perfect fluid approximation is 

= (E + P)uy.u v — Pg^u ■ (1.1.9) 

Using this expression and the definition of the energy-momentum tensor 

T IW = = P,x dpipd v ip - P<j, lv , (2.2.3) 


the energy E is obtained as E = 2 XP,x —P where the four-velocity u M 

To derive the adiabatic subtraction terms in following sections, it is need to obtain the 
expression of the sound speed. The sound speed of k-inflation is given by 




P, 


X 


,x 


E. 


2 XP,xx +P,x 


(2.2.4) 


can be negative and c s can exceed 1 by definition. However, c 2 s > 0 is required in general 
for stability [16] and c s does not exceed the speed of light as long as the Lagrangian does 
not contain the negative power of X. 

Also z is expressed in terms of P and X, 


c?\JE + P _ a 2 y/2XP, x 
c s T~L c s T-L 


2.2.2 The adiabatic subtraction term for the scalar perturbation 


The Mukhanov-Sasaki equation Eq. (1.1.21) is a second derivative equation and we can 
derive the adiabatic subtraction term for the comoving curvature (scalar) perturbation 1Z 
from this. 

( 1 . 1 . 21 ) 
( 2 . 2 . 6 ) 


v'k + ( c 2 s k 2 - =- ) v k = 0 


where 


v k = zTZ k , z = 


a 2 V2XP^ 
c s H ' 


cl = 


•X 


2 XR xx +P. 


;X 


In this case, using the notation in the section 2.1, 


2 2?2 

= c s k , <r = - 

z 

and the physical amplitude is schematically given by 


(|7?.(a:)| )phys — 


r°° rlk 

7^ k 2 [l^fc^lbare “ l^fcWlsub] 


where 


l^(’7)lsub = l^(^)| 2(0) + |A fc (r/)| 2(2) 


is a adiabatic solution of TZ up to second order. 

Then the adiabatic subtraction term for k-inflation model can be calculated 


M)\ 2 suh = ^ 


-1 -2 


1 _i d 2 -i 


1 


d?7 2 


+ 1(7 


1 


2 Z 2 Cek 


1 + 


2c? fc 2 


c?fc 2 


l _ c l 

4 c s 


3 c’ 2 
8 c? 


(2.2.7) 

( 2 . 2 . 8 ) 

(2.2.9) 


( 2 . 2 . 10 ) 


The third term in the parenthesis in the second line comes from the derivative of uj k . 
In the models having canonical kinetic terms (for example slow-roll model etc.), the sound 
speed c s is equal to one, and this derivatives are equal to zero. This result accords with the 
result in [3], which considered slow-roll inflation model. 
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2.2.3 The adiabatic subtraction term for the tensor perturbation 

The tensor perturbation or the gravitational wave hij is one of the inodes of the metric 
perturbations Eq. (1.1.7). The tensor perturbation is gauge invariant itself, and has only 
two degrees of freedom because of the property of the metric. 

While the scalar perturbation propagates at the sound speed, the tensor perturbation 
propagates at the light speed even in k-inflation case. The each degree of freedom of tensor 
perturbation obeys the following equation like Mukhanov-Sasaki equation with c s = 1. 

k 2 - v Tk = 0 (2.2.11) 

where vtu = f/ifc. h k denotes one of two polarized gravitational waves and we must add 
each results at the end of calculation. 

The subtraction terms can be derived similarly. In this cace, 

ujl = k 2 , <7 = -°- . (2.2.12) 

The calculation of the subtraction term has done by [3] in slow-roll model case, and the same 
analysis can be used in k-inflation model case. The expression of the subtraction term for 
tensor perturbation is obtained as 

M^lsub = ^ (i+ 2p—) - (2.2.13) 

and it is efficiently suppressed at late time. 


v Tk 


2.3 Non-minimally coupled k-inflation 


The inflation models which has the coupling terms of inflaton and scalar curvature have been 
considered. Higgs inflation model (see e.g. [18]) is one of these models, and of course there 
is no reason to prohibit inflaton from coupling to scalar curvature in general. Then let us 
generalize more the Lagrangian and obtain the adiabatic subtraction term for non-minimally 
coupled inflation models. 

The action of non-minimally coupled k-inflation model is given by 

S=\j d A Xs/—g [—f(ip)R + 2P(tp, X )] . (2.3.1) 


f(ip) = 1 corresponds to the minimally coupled case and f(ip) = 1 + \g> 2 corresponds to the 
conformally coupled case. 

In non-minimal case, the Einstein equation differs from the minimal couping case, and the 
scalar perturbation does not obey the Mukhanov-Sasaki equation. However, according to 
[19], it is shown that the scalar curvature obeys the equation like Mukhanov-Sasaki equation 
by using ADM formalism [11]. 


V 'k + ( C s,eff fc2 ~ — )vk= 0, V k = Z eS llk 

•2-eff / 


z e g and c S G ff have been estimated directly by ADM formalism in [19]. 


4r = 6e 2e 


n 




O' 

O' 2 - 9" 


- 1 +2 


. a 4 £ 

W 


- 9') 2 


a 4 £ 
p 2 6 


where 


9=-hi (fa 2 ), E = XP, x +2 X 2 P, xx . 


i.e., O' / 9 corresponds to P / H respectively in the minimal coupling case. 


(2.3.2) 


(2.3.3) 

(2.3.4) 


(2.3.5) 
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Using Eq. (2.3.2), we can get the subtraction terms for non-minimal coupled case, 


\K k (v)L h = 


2z. 


eff'-'S.e 


:ff k 


1 + 


1 


^eff 


1 


i a 


s,eff 


2cLfffc 2 Zeff ct_ Rff k 2 \ 4 c. 


!,eff 



(2.3.6) 


This subtraction term depends on time and /(</?) in a complex manner. Moreover, the 
ADM formalism uses a very long method of calculation to derive the equations of motion. 
Once we have the expression Eq. (2.3.6), there is no need to search for other methods. 
However, for checking physicality and simplifying the process, we will consider a different 
method to obtain the subtraction terms for non-minimal coupling case in Part 3 (based upon 
using the equations of motion). 


2.4 Time development of the subtraction terms 

2.4.1 Freezing out and the subtraction terms 

The bare power spectrum has a property whereby certain parameters “freeze out”: they do 
not develop after horizon crossing aH = k in canonical case, or the sound horizon crossing 
aH = c s k in the non-canonical case. The reason is that the comoving curvature perturbation 
1Z becomes constant in the large scale limit — ki] <1 or —c s kr) <C 1. 

The frozen modes of the perturbation become observable after the horizon re-entry, which 
occurs in radiation dominant era or matter dominant era, see Fig. (1). The arrows indicate a 



Scale factor 


Figure 1: The time development of the comoving Hubble radius 1 /aH (log-log graph). 

specific comoving scale 1/k of the power spectrum. The comoving horizon (or the comoving 
Hubble radius) shrinks during the inflation and the scale exits to the outside of the horizon at 
(a). After reheating around (b), the radiation dominant era begin and the comoving Hubble 
radius grows. The mode re-enter into the horizon at (c). It is considered that the primordial 
perturbations are transformed into the observables with some processes after (a). However, 
it is unclear when it occurred. 

The longer the time the perturbations are freezing, the more primordial information we 
can obtain. Therefore we are interested in the large scale limit. 

When the adiabatic regularization scheme was adopted to evaluate the power spectrum, 
it was considered that the subtraction terms also freeze out after crossing at first [1]. The 
subtraction terms for slow-roll inflation models at horizon crossing are not small and cannot 
be neglected, therefore the author claimed, see [1], that the amplitude of the perturbations 
was changed from the bare one. 

However, the time dependence of the subtraction term is not obvious from Eq. (2.2.10) 
even if we take the large scale limit. Moreover, it should be estimated up to the stage 
of reheating or up to times just before the primordial perturbations are transformed into 
classical quantities. 

The method of subtraction has been argued by various authors (see the review [9]). Based 
on the above, we adopt the scheme adopted by Urakawa and Starobinsky [3]. Their scheme 
is similar to Parker’s one, but not quite the same. In this scheme: 
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• We subtract up to the second adiabatic order term for the power spectra. 

• We consider the time development of the subtraction terms even after horizon crossing. 

• To consider time development, we do not require any cutoff of the range of the k 
integrals. 

The first condition is the minimal subtraction scheme as usual. The second condition is 
different from the original idea of Parker [1] , while the third condition is required to maintain 
the equation of motion of the inflaton in coordinate space [9]. 

Importantly, we are of the view that this scheme is simple and does not need any ficti¬ 
tious assumptions or processes; hence we use this method to analyze the observable power 
spectrum. 

In the following sections, we see how large a subtraction term remains after the inflation 
era. Although it may be not so important for observational physics, we also investigate the 
time development of the subtraction term during inflation in two specific cases, see Appendix 
B. 


2.4.2 The subtraction term in terms of the slow-roll parameters 

It is useful to rewrite the subtraction term in terms of the slow-roll parameters which we 
are familiar with to see the time development. There are several ways to define the slow-roll 
parameters. Among them, we use the Hubble flow functions and the sound flow functions. 
The Hubble flow functions is defined by 

_ din e n _ H[ n [ 

671+1 = dN ’ 

and the sound flow functions Si is defined by 

_ dlllS n ~ Cs,ini 

<Wl = j A r ) °0 ~ - 

dN c s 

where N = lna/ai n i is e-folding number. The subscript i„i means the initial value which is 
irrelevant, but we should fix it at the horizon exit to compare the subtraction term to the 
frozen bare spectrum. 

The definitions are also written in alternative forms 


(2.4.1) 

(2.4.2) 


£n+1 “ H e n ~ n e n ’ 71+1 “ H S n ~ U S n ' 

Then rewrite z in terms of <y. First, 

= i£° 

61 He o H 2 

This is equal to the slow-roll parameter e = defined in usual. 
From Friedmann equations Eq. (1.1.13) and Eq. (1.1.15) 5 , 

E = -3 H{E + P) = 6 HH 

E + P = -2 H 

Substitute this into Eq. (1.1.22), 


—2 H \/2e \cl 

= c s 


z = 


WE + P ay/E + P 


c s H 


JJ 


(2.4.3) 


(2.4.4) 


(2.4.5) 

(2.4.6) 


(2.4.7) 


Using this expression, we can rewrite z"f z [20]. 

— = a 2 H- /2 — e\ + —e 2 + ~ 2 £l£2 o 6263 + (3 — £i + £ 2)^1 + S 2 + (2.4.8) 


5 In first order perturbation theory, it is enough considering only the zeroth order part of 2 because the comoving 
scalar curvature 1Z itself is the first order of perturbation. 
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We can also rewrite the third them of Eq. (2.2.10). 


- — — - — —-a 2 H 2 5\(2 — 2ei + <5i + 2<5 2 ) (2.4.9) 

4c s 8 c, 8 

The calculation in detail is given Appendix A for completeness. 

Finally, we obtain the expression of the adiabatic subtraction term. 

^ fe(7?)ls2ub = 2^fc{ 1+ (S) ( 1 + <5e + ^)} (2A10) 

Here we defined two variables: 

6e = l + | e 2 - ^eie 2 + ^ 263 ) (2.4.11) 

Sc s = i<Ji (10 - 2 er + 4e 2 + 3<5i + 2S 2 ) (2.4.12) 

O 

and have isolated the explicit effect of the non-canonical kinetic terms. 

During inflation, we usually expect small slow-roll parameters, which can be constructed 
by the Hubble flow functions and the sound flow functions, and Se and Sc s are smaller than 
0(1). Even after inflation, Se is not so large because ei = 2 and e 2 = 0 in the radiation 
dominant era a(t) oc y/t, and ei =3/2 and e 2 = 0 in the matter dominant era a(t ) oc f 2 / 3 . 
How large Sc s is depends on the model. However, the “sound speed” c s does not exceed 1 in 
usual k-inflation models, and we expect that Sc s does not become too large. 

Therefore it depends on the behavior of the factors 2z 2 c k and ,-,~ 2 c k {jr^j how the 
subtraction term develops. Let us see how the each term behaves at late time. 

The zeroth order subtraction term 0 _ 2 C becomes exponentially small in general because 


1 c s 

2 z 2 c s k 4ei a 2 k 


(2.4.13) 


and a expands exponentially during inflation. e\ also become large while c s cannot become 
too large. This is a natural consequence because the zeroth order term corresponds to the 
zero-point energy in the Minkowski spacetime. 

Concerning the second order subtraction term, mainly the behavior depends on c s because 
it is known that (q^-) = 4 ^ fc3 at the late time is sufficiently suppressed compared to 
the its value at the horizon crossing in canonical case [3]. Because ei in post-inflation era is 
bigger than during inflation and H oc f -1 in radiation and matter dominant era. 

In the non-canonical case, 

1 faH 
2 z 2 c s k \c s k 

Therefore if there is a model in which c s becomes sufficiently small, the subtraction terms of 
the model can be large enough such that it affects the power spectrum. 


H 2 


4ei c s k 3 


(2.4.14) 
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Part III 

Conformal transformations and 
physicality 

In Part 2, we discussed how the adiabatic subtraction term for the non-minimal k-inflation 
model can be derived by using the ADM formalism. However, minimal coupling models 
are studied more often than non-minimal inflation models, where part of the reason is the 
calculations are easier. Therefore, it is convenient to analyze non-minimal models if we can 
apply knowledge from minimal inflation models in some form. 

Using a conformal transformation, we can rewrite the Lagrangian of the non-minimal 
coupling model as that of one in the minimal coupling model. The metric and the pressure are 
altered by such a conformal transformation, hence the physical properties of the conformal- 
transformed universe are different from the true minimal coupling case, even though the 
action has the same form as a minimal coupling model. 

The frame which has the non-minimal coupling term is denoted the Jordan, while the 
conformally-transformed frame is denoted the Einstein frame. It is interesting to note that 
the physical interpretation is completely altered, yet the observations do not change in the 
transformation (see e.g. [21]). 

Importantly, in this chapter, we will show that we can derive the adiabatic subtraction 
terms (extracted from the Mukhanov-Sasaki equation) for non-minimal k-inflation with the 
conformal transformation, without using the ADM formalism. 


3.1 Einstein equation for non-minimal coupling model 
in Jordan frame 

Before the discussion of the conformal transformation, let us review the calculations in Jordan 
frame. How different it is from the minimal case is also argued. 

The action in the Jordan frame is 

8=\J d*x^g[-f(<p)R + 2P('p, X)] (2.3.1) 

and the Einstein equation becomes 

5 ^<>VaV a / - V M V„/ + / ^Rfiv - * Rg^ - P,x + Pg ^ = 0 • (3.1.1) 

To simplify, define three tensors and rewrite the Einstein equation as follows: 


G fils - RflU 2^t JLV 

(3.1.2) 

Tfiu = P->X ^ 

(3.1.3) 

Fnv = 3^i/VaV a / - V M V^/ 

(3.1.4) 

fG» v = n - n 

(3.1.5) 

If we pick up only the first order part of the perturbations, 

S(fG^) = Gfjf + fcSG^ = ST% - 6 Ft . 

(3.1.6) 


Hereafter the subscript c denotes the classical part and 6 denotes the first order part of the 
perturbations respectively. 

The gauge fixing is needed to calculate the equation of the scalar perturbation. We use 
the conformal-Newtonian gauge Eq. (1.1.17) in this section because the correspondence to 
gauge-invariant variables is clearer than other gauges. 
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From the transformation law of the tensor, the gauge-invariant <5G M „, ST ll v , and 5F^ V is 
written by 


<5G° 0 = 6G° 0 - (G c ° o y(B - E') 

SG'j = 6G i j -(G c i j )'(B-E') (3.1.7) 

) (B - E’) 

and so on. The gauge invariant Sf is also obtained from the transformation law of the scalar, 

Sf = Sf- f' c (B - E) . (3.1.8) 


SG , = <5G ,■ - ( G c ° 0 - -G 


In conformal-Newtonian gauge, B = E = 0, and we get the gauge-invariant Einstein 
equation by replacing the tensors by the gauge-invariant ones. 

g c v 7+/ c *n=- sr v ■ (3-i-9) 

Then calculate the tensors and the equation in accordance with the definitions Eq. 
(3.1.2)-(3.1.5). After the tedious calculations, we can obtain SG [22] 


5G° 0 = 

4{A^-3^(vE' , + ^$)} 

a z 

(3.1.10) 

II 

O 

4(^'+H$),i 

or 

(3.1.11) 

<• = 

2 

jV' + U{29 + $)' + (2 U' +U 2 )<& + ^ | $ij 



i-1 

ST 


(3.1.12) 

where A = and ST 41 v . 



6T°o 

= 2 XSP, x -SP 

(3.1.13) 

ST ° 

= ~i;Pc,xv'S<Pn 
a z 

(3.1.14) 


= -SPS% 

(3.1.15) 


The Einstein equations are then as follows: 

+ = ^a 2 (5T° 0 ~5F 0 0 ) (3.1.16) 

fcV&' + HQU = \a\HT\-JF\) (3.1.17) 


(n 1 + ^n 2 j sfSij + / c j^ ,, + ^(2^ + $) , + (2^ , + ^ 2 )$+^A($-*)J< 5 iJ • 

- \fc(* - *),« = - SF 1 ,) (3-1-18) 

Let us set i ^ j in the Eq. (3.1.18). Because the non-diagonal content of ST',- is zero by 
definition then 

f c (*-*),ij=-a*6F i j . (3.1.19) 

The spatial non-diagonal components of the Einstein equations are zero in the minimal 
coupling case. The existence of this anisotropic inertia is a big difference between the minimal 
and non-minimal couplings model. 

The explicit form of depends on the model, but we calculate the general form for 
generality. 
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From Eq. (3.1.4), 

n = Kgf (d p d x - ( 0 ) r« A a a ) f c - g pp (d p d v -^v%d a ) f c 

+ n ( g Xp d p d x Sf + 6g Xp dpd\f c - g xp ^T a pX d a Sf - g Xp ST a pX d a f c - 5g Xp ^T a pX d a f c ) 
- grdpdjf - 5g pp d p d„f c + g ppW T a pu d a 6f + g pp 5T a p „d a f c + Sg pp ^T a piy d a f c (3.1.20) 


where the Christoffel symbol T“ A =^°7“ A + JT“ A (the notation is temporarily changed for 
simplicity). The first line is the classical part F r , and the rest is the linear order part SF 


In up to first order perturbation theory, SF p v (5f) = 

SF l f(Sf) and we get 


5F °0 

= 4 

a z 

(-A Sf + 3HSf' - 6 H<S>f' c - 

3*7') 

(3.1.21) 

SF°i 

= 4 

a z 

(-Sf + HSf + Qf^ ti 


(3.1.22) 


1 

s\ ? {sf - A Sf - 2 4>/" + nsf - (4> + 24/)'/' - 2 W'} 



~\~d 2 



(3.1.23) 

Then when 

* 7 i j, 






/ c ($- $),y = -a 2 IF\ 

= -didjSf 

(3.1.24) 


In the minimal coupling case, there is no anisotropic inertia and we can combine the Ein¬ 
stein equations to get the Mukhanov-Sasaki equation [22]. However, due to the anisotropic 
inertia term and other additional terms from Sf, it is difficult to combine straightforwardly 
the equations in the non-minimal coupling case without resorting to different mathematical 
techniques, such as the ADM formalism. 


3.2 Conformal transformation and invariance 


3.2.1 Conformal transformation 

Again consider the non-minimal coupling k-inflation model in the Jordan frame. 

S = \J d 4 xV^g [~f{<p)R + 2 Pfr, X)] (2.3.1) 


To derive the adiabatic subtraction term without ADM formalism, perform the conformal 
transformation g pv —> fj pl , = f(y>) gp , v to obtain the action in the Einstein frame [23, 24]: 


where 


S = 


1 

2 


-R + 2P(ip,X) 


a(jj) 

R 

X 

P 


- f R-3V p Wlnf-*(d p lnf)(d p lnf) 
^g pv d p yd v y = jX 
^P+^- f (d p \nf)(d p \nf) . 


(3.2.1) 


(3.2.2) 


Note that the coordinates are not changed by the conformal transformation when we use 
conformal time, i.e.) drj = di 7 , dx l = dx l . Then the line element in the conformal-Newtonian 
gauge becomes 


ds 2 = f(ip)ds 2 

= f((p)a 2 (i 7 ) [(1 + 2^)dg 2 — (1 — 24')c&? 2 ] 

= d 2 {rj) [(l + 2 $)dr ? 2 -(l- 2 $)<ir 21 


(3.2.3) 
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After expanding each coefficient up to linear order of perturbations, we can get the transfor¬ 
mation laws of the scalar perturbations [25]. 

If ~ If , 

$ = $ + —, T = ^-- (3.2.4) 

2/c 2 f c V ' 

Let us check whether or not the anisotropic inertia in the Jordan frame vanishes in the 
Einstein frame. Substituting Eq. (3.2.4) into Eq. (3.1.24) leads to 


fc 




,ij = fc - $) nj -didjdf 


—didjSf . 


(3.2.5) 


Therefore, we can find the correct relation: 4> = 4/ in the Einstein frame. 


3.2.2 Conformal invariance and physicality 

The conformal invariance of the comoving curvature perturbation 1Z and its correlation func¬ 
tions (7 Z ...) was shown in [24, 25]. i.e., 


H = k, (n...) = (iz...) 


(3.2.6) 


Then, if we regularize the correlation functions by adiabatic regularization, the subtraction 
terms in each frame must also agree with each other. Because otherwise the regularized 
power spectra depends on the frames and the method of either conformal transformation or 
the adiabatic regularization is spoiled. 

In this section, we will derive the adiabatic subtraction terms in the Einstein frame and 
show the conformal invariance of the adiabatic subtraction terms in both the Jordan and 
Einstein frames. 

First, consider the equation from the action Eq. (3.2.1) 


d 2 ^ 2 2 1 d 2 z\^ 

—v k + c 2 s k 2 - ) v k = 0 

dr/ z \ z dij z / 


where 


Vk 


= zR k =zlZ 


k, z = 






’X 


2XP,2x +P>x 

From the universality of the conformal time, the equation becomes 


Vk 


c s k 2 - — ) v k = 0 . 


(3.2.7) 


(3.2.8) 


(3.2.9) 


We can get the adiabatic subtraction terms for the power spectrum from Eq. (3.2.9) as the 
same way in the minimal coupling case. 


\nv)\i h = 


i 


2 z 2 c,fc 


1 


+ 


1 


2 zcik 2 elk 2 


1J 

4 c„ 


3? 


^ 2 


(3.2.10) 


From the adiabatic subtraction term in the Jordan frame Eq. (2.3.6) and in the Einstein 
frame Eq. (3.2.10), it is shown that the adiabatic subtraction terms of non-minimally coupled 
k-inflation model depend on only z e g / I and c Sie ff / c s . Therefore, if z and c s are equal to 
the effective quantities Eq. (2.3.3)-(2.3.4) after conformal transformation, the subtraction 
terms in both the Jordan frame and Einstein frame correspond to each other. 

Once it has shown the correspondence of the subtraction terms in both frame, it is not 
necessary to derive it in Jordan frame by using ADM formalism. We can calculate it in 
the Einstein frame, and transform it into the Jordan frame if we want to see the physical 
property due to the non-minimal coupling. 

For conformal transformation, it is convenient to define two functions like Eq. (2.3.5). 

0 = \ hi (fa 2 ) = i In (a 2 ), E = WPy + 2X 2 Pyy (3.2.11) 
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Using Friedmann equations in the Einstein frame, c s and z are expressed by 


XP,x 


V II 2 


z 2 = 


2 ei a 2 


E E 

Thus we should transform three quantities H , el and E. 


2 a 2 E 
H 2 


if = 


ei = 


E = 


S = ^=e-V 

a a 2 

= 1 _ # = 1 _ 0" 


0 

/ v 0 * 


H 2 


9' 2 

3X /dln/\ 2 ' 
2 f \ dip ) 


C2 ^ ^ 

IV / . 5 \ 2 f 1 „ 3X fd\nf\ 


+2( /J vW 1 r- P+ 2 / l V J 


J 2 XP ’ x +2X2p ’ x * + 2 fX ) 


(dlnfY 


«“i(f 


= e _ 40 a 4 


o 20 


= e 


—40/V2 


£ + 3^(0' 


3e- ( P - 1 


a 4 E 

6»' 2 


To transform E, we used Eq. (3.2.2) and the following relation. 


P = L P+ ^( dll lnf)(d v \nf) 

■ y-MW 


^(^ln/X^ln/) 


Substitute them into (3.2.12), 


= 2 i 3e 28 I U ~ 1 


„20 


b 2 = 


9' 


1 _ 

-L n / 2 




a 4 E 

~9' 2 ~ 


9' 2 - 9" 


a 4 E 


3 e 2 «(«- l ) 2 + ^ 3 ( H - 0') 2 + $ 


(3.2.12) 


(3.2.13) 

(3.2.14) 


(3.2.15) 


(3.2.16) 


(3.2.17) 


(3.2.18) 

(3.2.19) 


Comparing the results with Eq. (2.3.3)-(2.3.4), it is realized that the effective variables in 
the Jordan frame and the variables in the Einstein frame agree with each other. 


— Z , C Sj eff — C s (3.2.20) 

Therefore the adiabatic subtraction terms in both Einstein and Jordan frame are identical, 
and we get the same results of regularized power spectrum in either frame. 

Of course, because the comoving curvature perturbation itself is conformally invariant, 
the equations the perturbations obey should also be conformally invariant. The above cal¬ 
culations show this explicitly. 
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Part IV 

Conclusion 


We have investigated the adiabatic regularization of the power spectrum in general single 
scalar field inflation. We derived the adiabatic subtraction terms and found that (due to the 
sound speed) their behavior is not trivial. The analysis was partly done by assuming some 
specific models or approximations [26] , and there may be the models which have growing or 
non-vanishing subtraction terms. 

We have also calculated the subtraction terms for the non-minimally coupled inflaton 
case. The subtraction term is completely determined by the equation of motion of the frame- 
invariant scalar curvature perturbation, and we confirmed its frame invariance explicitly. 


Discussion and future work 

What do vanishing subtraction terms mean? 

In section 2.4, we have considered the time development of the adiabatic subtraction terms. 
We saw that the speed of sound may be able to make the subtraction terms large enough 
to affect the observable power spectrum: in canonical inflation, it generally decays after 
inflation era (see Appendix C) and we expect that it becomes sufficiently small when they 
are transferred to observables. 

We assumed a scheme where the bare spectrum freezes at the (sound) horizon crossing, 
while the subtraction terms remain time-developing. The vanishing subtraction terms mean 
that the integral of the bare spectrum at the crossing has no divergence to begin with. This 
result agrees with the conventional interpretation about the observable power spectrum: the 
observable spectrum has no high frequency modes because they have not yet exited the 
horizon and not yet become “classical”. 

However, even if the result with regularization does not differ from the result “without” 
regularization, it does not mean that the argument for regularization of the power spectrum 
is unnecessary. Indeed, it gives mathematical support to the observable power spectrum, and 
also it is needed in an interacting theory. 

Future work 

In section 2.4, we have investigated how the subtraction terms behave qualitatively at late 
times with the subtraction scheme used in [3]. Related to this scheme, there are some 
problems still to be solved. 

The subtraction scheme was not established completely. There are many discussions 
about this issue (e.g. [3, 9, 27-29]) and some papers such as [9, 27] claim that the regular¬ 
ization of the power spectrum is originally unnecessary. Some of these claims have already 
been rejected, nevertheless, the argument is still continuing. This is possibly due to our poor 
knowledge of the “freezing out” and the unclear nature of the renormalization condition 
using adiabatic subtraction (see also [30]). 

Understanding “freezing out” is an important problem of cosmology. While it is true that 
the comoving curvature perturbation loses its time-dependence at large scales: — krj <C 1, 
how to become classical is still unclear. The scheme we take in this thesis assumes that 
the freezing out does not occur at the horizon crossing, in a sense. When the quantum 
fluctuations become observables is a fundamental problem we should consider in more detail. 

We should also mention that we have only calculated the value of the subtraction terms 
qualitatively. As we saw in section 2.4, it is possible for terms to survive or become large 
when the model has non-canonical kinetic terms. We therefore need to know how large 
these surviving terms are and compare the regularized power spectrum with the bare one. 
Some models might need to be modified to achieve a realistic power spectrum if the sub¬ 
traction terms remain significant at the time of classicalization: we can use the regularized 
power spectrum to constrain an inflationary model by combining it with other theoretical 
constraints. 
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Part V 

Appendices 


A Hubble flow equations in detail 


The effective potential of Mukhanov—Sasaki equation 

Let us derive the explicit form of the effective potential Eq. (2.4.8) 

z" ( 3 111 

— = a 2 H~ < 2 — ei + —£2 + 7£2 ~ 76162 + 76263 + (3 — ei + £ 2)^1 + $1 + < 5 i <$2 

Z A 4 A 

From Eq. (2.4.7), 
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Then let us calculate the time derivative of z. After some calculation, we get 
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Substitute them into Eq. (A.l). 
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Then rewrite this in terms of the Hubble and sound flow functions. Eq. (A.4) contains 
the third derivative of H and we need up to the third Hubble flow function £ 3 . 

The explicit forms of the flow functions are below. 
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2 _fr 2 A can rewrote by using £ 263 , but then the term proportional to H2 ^ 2 
remains. To rewrite it we need the square of £ 2 - 
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Rewrite the second line of Eq. (A. 4) . 
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Also up to the second sound flow function 82 is needed. 
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Rewrite the third line of Eq. (A. 4). 
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By combining Eq. (A. 10) and eq.(A.13), we can get the explicit form of the effective 
potential Eq. (2.4.8). 


The third term of the adiabatic subtraction term 

Next, let us rewrite the third term of the subtraction term Eq. (2.4.9). 
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To use the above calculation, we need translate this terms in terms of proper time deriva¬ 
tives. 
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Combining Eq. (A. 14) with the effect by time dependent sound speed in the effective 
potential Eq. (A. 13), we obtain 

— {(3 — t\ +£ 2)^1 + 8 \ + ^ 1 ^ 2 } — (2 — 2ei + <5i + 262 ) = — <5i(10 — 2ei + 4e 2 + 3<5i + 282 ) • 

2 00 

(A.15) 

Then we have rewritten the adiabatic subtraction term in terms of the flow functions 
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B The time development of the adiabatic subtraction 
terms during inflation 

We have investigated the time dependence of the adiabatic subtraction term for non-canonical 
inflation in the section 2.4. The zeroth adiabatic order term decays in general, while we 
cannot say anything about the second order term unless we specify the models (see [26]). 

From the view point we take in the section 2.4, the time dependence of the subtraction 
term during inflation is irrelevant. However, some people still claim that the subtraction 
term should be estimated at the horizon exit [15, 31]. 

Hence we estimate the time development of the subtraction term during inflation in the 
following two cases: (1) the slow-roll inflation model in quasi-de Sitter spacetime and (2) 
DBI inflation model. As a result, we find that the subtraction terms at the horizon exit 
exceedingly suppress the bare power spectrum in these two cases. 


B.l The slow-roll inflation model in quasi-de Sitter spacetime 

Let us review the case of slow-roll inflation model in de Sitter background spacetime first. 
Of course, however, the exact de Sitter universe can be achieved by the cosmological con¬ 
stant rather than inflaton. We set the slow-roll inflaton with perturbations on a de Sitter 
background, and hence the derivation of the power spectrum is an approximation. 

In the case of exact de Sitter universe, 

a(t) = e Ht — a(ri) = (B.l) 

and H is constant. The Hubble flow functions e» (i > 1) are equal to zero by definition. The 
sound flow functions are also zero. Furthermore, we are assuming a model with canonical 
kinetic term. 

Then the effective potential of Mukhanov-Sasaki equation becomes 


— = 2 a 2 H 2 = 4 (B.2) 

z rf 

and the following expression is obtained. 

(l + ^2) = ^3 ( fc V + !) (B- 3 ) 

Here we use that E + P = p 1 . In the large scale limit, 

_Km o |K t (,)|, 2 „ b = (B.4) 

This is equal to the bare power spectrum of the slow-roll model in the large scale limit. 

Under the slow-roll condition, we assume that the second derivative of inflaton p is small. 
We can neglect the time dependence of p in the denominator of Eq. (B.3) during the inflation, 
and find that the first term exponentially decays but the second term is constant in the time 
development. It is shown in the Fig. (2). 

Take the difference between the bare power spectrum after the horizon exit (in large scale 
limit) and the subtraction term, 
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(B.5) 


Because H and p is nearly constant during slow-roll inflation, the physical power spectrum 
in large scale limit is exceedingly suppressed by the subtraction term. This result agrees 
with the prior research [1, 15, 31] which argue the amplitude of the fluctuation of not the 
comoving curvature but the inflaton itself: (| <5 v? | 2 ) P h y s is also suppressed by the subtraction 
term during inflation. 
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Figure 2: The time development of the subtraction term Eq. (B.3) in the case of de Sitter spacetime. 
We are assuming <p = 1 (constant) to plot. We also set H = k = 1 because these constants affect only 
the overall magnitude of the power spectrum. The subtraction term does not become zero even at the 

late time of inflation era. 


B.2 DBI inflation model 


Next let us to see how the adiabatic subtraction term behaves in DBI inflation model [32]. 
DBI inflation is motivated by string theory and it achieves the inflation by the “D-cceleration” 
mechanism. The cosmological property of this model has been analyzed, see e.g. [32, 33]. 
The Lagrangian of DBI inflation model is given by 


P{y,X) = (y/l - 2M<P)X - l) - V{ip) 


(B.6) 


where /d(^) is the (squared) warp factor and we redifine X in terms of the proper time 
metric 6 , i.e., X = \<p 2 . 

The adiabatic subtraction term in this model can be calculated by using the general result 
Eq. (2.2.10) or Eq. (2.4.10). 
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Let us calculate the speed of sound and z. 
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Note that both c s and z do not depend on the potential V ( ip ) in general (not limited to DBI 
model) because E + P is derivative of X'. However, the behavior of tp is determined by the 
f D M and VM- 


b We set the sring coupling g s = 1 because it only affects the overall magnitude of the power spectrum at large 
scale limit. 
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Figure 3: The time development of the subtraction term in the case of a DBI inflation model fo{v) = A) 
V{ip) = m 2 <p 2 . We set k = 1, ejj = 0.01 and m = 10 -15 to plot. The subtraction term decreases very 
rapidly but the second order term which cancels the bare spectrum remains. 


We can obtain the subtraction term in terms of /d(v?) and X by substituting 2 and c s 
into Eq. (2.4.10). However, there are many choices of the two functions. In this section, 
we assume a simple model which has foif) = A/<p 4 and V(ip) = m 2 <p 2 . A is the ’t Hooft 
coupling and m is the mass of the inflaton. This model is analyzed in detail in [33] and can 
achieve a power law inflation. 

In this model, each function behaves at late time (but it is still inflation era) as following 

[33]. 
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where ejj is a slow-roll parameter. 
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i.e. c s <C 1 is necessary to achieve accelerated expansion. 

e d can be calculated by using the Friedmann equation and Eq. (B.ll) 


bo 


3 

\J 3m' 2 A + 2>/3m 2 A 



(B.13) 


and it is approximately constant at late time. Let us see the relationship between eo and 
the flow functions. 
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Therefore eo = e\ = — ^ (5-| and Se and Sc s are constants during the inflation. 

Then we can also use the factors 9 ^ c k an d 0 , 2 C k (fr) to see ^ 3me dependence 
of the subtraction term. Substitute Eq. (B.ll) into them and we obtain the following 
expression. 
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The zeroth adiabatic order factor becomes small, while the second term is constant. This is 
also the same as the bare power spectrum of this model is given by (in the Bunch- Davies 
vacuum and in large scale limit) 

1 TJ1? 

I^fc^lbare = ^4^3 = 6^3 ' ( R18 ) 

The time development is shown in the Fig. (3). 

Therefore the regularized power spectrum at the large scale is also strongly suppressed 
in this model. 
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C The time development of the adiabatic subtraction 
terms after inflation 


As mentioned, we assume that the subtraction terms do not freeze out at horizon crossing. 
The subtraction terms shrink after inflation in the canonical kinetic term case [3]. In the 
non-canonical case, they also decrease if we assume the nearly scale invariant power spectrum 
[26]. 


Then, the following question comes to mind; How small the subtraction terms become? 
No one seems to check this quantitatively. In this appendix, we see it via a very rough 
estimate for the following research. 

We restrict ourselves to considering the models which have neither the non-canonical 
kinetic terms nor the coupling with the scalar curvature. Also we assume that the scalar 
perturbation obeys the only one equation throughout the history of the universe, and we 
neglect the dark energy dominant era. 

The second adiabatic order subtraction term is 
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H 2 
4ei fc 3 
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For the simplicity of the formula, we set the magnitude of the dimensionless bare power 
spectrum to 1/87T 2 . Then, to compare the magnitude of the subtraction terms with the bare 
one, we only have to see the time development of Eq. (C.l) for a specific mode. Eq. (C.l) 
can be rewritten in terms of slow-roll parameters [3]. 


H 2 _ i 
4ei/c 3 4fc 3 6 ^ 
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dV(2ei(iV) + e 2 (V)) 
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where N * is the e-folding at the horizon exit. Here we define the e-folding N(a) = In (a/a*) 
with ao = 1. 

We need separate the range of integral to calculate this. We consider three parts: inflation 
era (from N* to N enc j), radiation dominant era (from _/V e nd to N eq ), and after that. 

We assume the inflation ends abruptly at N en d. The behavior of the subtraction terms 
during inflation era can be calculated (see Appendix B). It has model dependence through 
the slow-roll parameters and e-folding, but is typically reduced to 10 _1 -10“ 3 . 

After the inflation, the slow-roll parameters become large and the time behavior of the 
subtraction terms lose the inflational-model dependence (except for the sound speed). 

In radiation dominant era, the integrand of Eq. (C.2) is equal to 4. The subtraction 
terms become small exp (—■4(iV eq — N en d)). N eq — N end is depend on the temperature at the 
end of inflation era. However, using the scale factor at the recombination and at the time 
the energy density of radiation and matter became equal, the subtraction terms are reduced 
at least to 10 -2 from the value at the end of inflation. Because 


N eq - N end > In (« 0.8 (C.3) 

\ UCMB / 

where ocmb is the scale factor at the recombination. 

The ambiguity of the calculation becomes trivial in matter dominant era. In this era, the 
integrand of Eq. (C.2) is equal to 3 and 

N 0 -N eq = ln(—) «6 . (C.4) 

\a eq J 

the subtraction terms are reduced to about 10~ 8 . 

Combining these results, the adiabatic subtraction terms are reduced at least to 10 -11 by 
the present. However, the power spectrum including subtraction terms should be estimated 
at the time it become “classical”. 

It is unknown when and how the power spectrum became classical. It perhaps is in 
inflation era or radiation dominant era. If the subtraction terms was not sufficiently small 
at the time, we may have the opportunity to observe the difference. 

This calculation is very roughly. We need more precise consideration including the effect 
by reheating and higher order quantum corrections. 
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D The subtraction terms in the Jordan frame 


In this section, we calculate the explicit form of the subtraction terms in the Jordan frame 
by using the conformal transformation to see how the non-minimal coupling term affects on 
it. 

The subtraction term has been derived in the section 3.2.2. 
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Let us rewrite z in terms of z which is the variables in the minimal case. 
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To derive the explicit form of the subtraction terms in the Jordan frame, we need the 
second derivative of z and c s . We can use some convenient calculations as below. 
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Derivatives of £ 

Let us calculate the derivatives of z first. For simplicity, define a function a, 
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(D.4) 


This function approaches to zero in minimal coupling limit bacause 9' is equal to B in the 


case of minimal coupling model. 
Then z 2 becomes 
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Again, to simplify this, define /? vanishing in minimal coupling limit 


„ , ( 2 2 v|«' 6" 

P = ota +(z + a ) ( — — 


so that 

Similarly, 

d 


z ,2 = 


1 / H 


(j) izz'+P? ■ 


(D. 8 ) 


(D.9) 


drj 


~,2\ 2 /HV, , 2 \(W 0 "\ zz' + aa' , ,z ,2 +^"+£' 


2 2 + a 2 \9 


% 8 ') z 2 + a 2 zz' + P 


(D.10) 


-y, a 5*2 2 rw 

dr? 

T “ J tZ 2 

dr] 


U(«' + s 2 {(¥ 


j _ ■Z 2 : / +Q'Q: / _i_ 

z 2 +a 2 ' -Z2/+/3 


: } 


2 © {(* 2 + « 2 )(¥-?)+«'+»«■} 


zz' +P 

z 2 + a 2 


(¥ 


;;'+Qa / i z ,2 +zz"+/3' 
0' ) " ; 2 +a 2 " l " zz'+,3 


2 W 0" \ , zz'+aa' 


(D.ll) 


n e> J ^ z 2 +a 2 


This is the expression of coefficient function of the second term of Eq. (3.2.10) in the 
Jordan frame. 


Derivatives of c s 

Next, we need the expression of the derivatives of c s in the Jordan frame. By definition 
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Using Eq.(D.5) and Eq. (3.2.14), 
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and it vanishes in the minimal coupling limit. 1 — ^2 — ei is the first Hubble flow function 
of the minimal coupling model. 

Then we get 
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and calculate the derivatives. 
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Define k as 
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where ft is defined as Eq. (D.18). 
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Results 


Substitute the expressions to Eq.(3.2.10), we get 
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We have to calculate the derivatives of a, /3, and k to obtain the final expression. How¬ 
ever, the expression is very complicated and there seems to be no further use in continuing in 
this fashion. It is actually enough to see that the coincidence of the expression in the minimal 
coupling limit, and the results definitely coincide. Although our expression of the adiabatic 
subtraction terms for “general” (non-minimal) k-inflation depends on the non-minimal cou¬ 
pling term in a complex manner, it includes the expression for the usual (minimal) k-inflation 
as a special case. 
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